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ABSTRACT 



It is observed that normally, during the polar night, 24- 
hour periods with a clear sky and light to moderate winds have 
their lowest temperatures in the middle of the period, and the 
highest temperatures occur at the beginning and end. For 
cloudy and windy periods, the shape of the temperature curve 
is reversed, with the maximum temperature occurring near the 
middle of the period. Despite many efforts by various invest- 
igators no physical explanation could be found for this vari- 
ation. 

Hisdal discovered that for clear skies and low wind speed, 
the same shape of temperature variation curve occurs, whether 
the starting time is 0000 hr, 1200 hr, or 1400 hr. Thus, the 
abnormal behavior of the temperature variation appears to be 
influenced primarily by the statistical bias involved in the 
selection criteria. 

This investigation, using six years of surface data, di- 
vides the diurnal temperature variation for various special 
types of days into two parts, the meteorological or diurnal 
variation and the superimposed statistical bias. Finally, the 
meteorological diurnal temperature variation, based upon all 
observations during the six year period, is determined. 

The writer wishes to express his appreciation for the 
advice and guidance given him by Professor W. van der Bijl of 
the U. S. Naval Postgraduate School in this investigation. 
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Introduction 



In his meteorological investigations with the Norwegian- 
British - Swedish Antarctic expedition of 11949 ■= 1952, V, 
Hisdal [l ] discusses a much pondered peculiarity of the diur- 
nal temperature variation during the polar night. For the 
various polar stations investigated, days with clear skies 
and/or light to moderate winds had their mlnlmiim temperatures 
near 1200 GOT, whereas for the cases of overcast skies and/or 
windy days the opposite type of diurnal variation occurred; 
that is, a maximum temperature near 1200 GOT. 

Concerning this strange temperature behavior, Simpson 
[2, p.7^ wrote in 1919s 

I have spent more time over this problem than over 
any other discussed in this work the investi- 

gation has led only to negative results. 

In 1951 Court [3] could still not explain this tempera- 
ture behavior. He wrote? 

Why should the maximum temperature on clear days 
during the antarctic night occur, almost without 
exception, in the hour after midnight? Simpson 
pondered this problem longer than any other of 
the many he investigated in compiling the meteor- 
ological reports of the last Scott expedition, 
without developing a satisfactory hypothesis. 

Rouch l 4J also found a daily temperature variation during 

the dark season, based on an analysis of the data both from 

the Charcot expeditions and from the Bird expeditions. He 

wrote that the average daily variation 

is a simple curve with a maximum around 1 or 2 
A.M. and minimum after noon .... I carmot give 
a satisfactory explanation. 
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Hlsdal's explanation for this phenomenon was simply that 
a physical explanation was not necessary at all. It could be 
explained in terms of an "automatic'' daily variation which re- 
sults from the nature of the type of day selected. For ex- 
ample, Hisdal [5, p.104] has shown that if a type of day is 
characterized by the presence of a temperature inversion, it 
will usually have a lower average temperature than the non- 
inversion days. This means that if the average diurnal tem- 
perature curve is computed using only the inversion- type days 
(having lower temperatures), the non-inversion days (having 
higher temperatures) will automatically tend to give a contri- 
bution near both edges of the 24-hour time span. The result 
of this automatic dally variation, shown schematically in fig- 
ure 1 , is that it superimposes a statistical bias upon any 
given diurnal variation that might be present. Moreover, the 
magnitude of the bias could be large enough to mask the diurnal 
variation completely. 

An interesting fact is that Hisdal 's explanation for the 
peculiar behavior of the diurnal temperature variation is ap- 
parently an independent discovery of the theory that Bartels 
C6H described in his Gottingen dissertation of 1922. 

Hisdal demonstrated this statistical bias explicitly by 
considering two diurnal temperature variation curves for 
polar stations, based on data for dark season days with wind 
speeds less than 15 miles per hour and cloudiness less than 
five- tenths. The first curve begins at 0000 GOT and has a 
minimum temperature at 1200 GOT. The second begins at 
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Time 

(b) 




Figure 1 

Schematic diagram showing the statistical bias effect 
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1200 GOT and has a minimum at approximately 0000 GOT. If 
there were a uniform diurnal temperature variation, the min- 
imum temperatures should always occur at the same local time. 
Moreover, varying the Initial time of the graph should not 
change the time of occurrence of the minimum. A logical ex- 
planation of the difference In these diurnal temperature vari- 
ations Is that the minima were primarily due to the statistical 
bias, not by dlumally different physical processes. 

W. van der Bljl [7, p.l49l also supports the theory of a 
"statistical bias" In his investigation of the diurnal pressure 
variation In northwest Europe. His study discloses a similar 
automatic diurnal variation due to data-selectlon bias. The 
basis of selection was that each day had to be In the month of 
December and have at least 5.7 hours of Insolation. These cri- 
teria resulted in a diurnal trend with low pressure at the 
beginning and end of the period, with high pressure In the 
middle. 

Thus, It has been shown that the existence of a "statis- 
tical bias" or "automatic diurnal variation" gives a reason- 
able explanation of the behavior of the temperature trend at 
polar stations during the dark season. However, Hlsdal In 
his investigation did not deduce the mean meteorological diur- 
nal variation of temperature, even though a knowledge of the 
effect of the statistical bias made this deduction possible. 

The purpose of this study is to compute the mean meteor- 
ological or diurnal temperature variation during the polar 
night. As a secondary purpose, the temperature variations 
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Induced by the statistical case-selection procedure will also 
be computed. 
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Data and Processing 



The data used in this investigation were obtained from 
the National Weather Records Center, Asheville, North Carolina. 
The meteorological information consisted of punched IBM cards 
of surface data at three-hourly Intervals from McMurdo Sound, 
Antarctica (78S, 166 b), for the years 1956 through 1961. Fig- 
ure 2 shov/s the geographical location of the area concerned. 

The duration of the polar night at McMurdo sound was de- 
termined by the use of the solar altitude curves prepared by 
LeVasteau LSj . These indicated that the duration of the dark 
season was from 20 March through 20 September. However, this 
period does not eliminate twilight effects and incoming solar 
radiation due to unusual atmospheric refraction conditions. 
Therefore, the shorter period, 15 May through 15 August, was 
used in order to Insure complete darkness. 

In processing the punched data cards, the IBM 82 card 
sorter was used to select observations for the 15 May through 
15 August period. Then the card data was transferred to a 
magnetic tape by use of the IBM 1401 printer. The magnetic 
tape was subsequently used in conjunction with a program for 
the Control Data Corporation 1604 digital computer, in order 
to determine: 

(1) the mean diurnal variations of temper- 
ature, cloudiness, atmospheric pressure, 
and wind direction for various years; 

(2) the same as (1), except for various 
starting times of the period; 
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(3) the same as (2), except for periods 
other than 24 hours; 

(4) the total six-yearly mesoi variations 
as stated in (1), (2), and (3) above; 

(5) and the Lament correction for the above 
mentioned curves. 




Figure 2 

Geographical location of McMurdo Sound, Antarctica 
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The Lamont correction 



The Lamont correction was mentioned In section 2 as one 
of the computer subroutines. The purpose of this correction 
Is to make the starting value of a continuous parameter and 
the final value of the parameter equal to each other. This 
technique simplifies the study of any cyclical variations and 
Is especially valuable In determining a diurnal variation. 

A more complete discussion of the Lamont correction is 
presented by Conrad and Poliak , For this study the cor- 
rected temperature was computed, using equation (2): 



"Hl - Tj 


(2) 


q; = l(To-Tj 
— _ — 


(3) 



where 



T;l = corrected temperature 

T: = sequence of three-hourly temperature observations 
(i — 0, 1, 2, •«•••«. .s.o.s, n) 



gj = Lamont correction. 

A typical application of the Lamont correction follows. 
For this purpose, use Is made of equations (2) and (3), and 
the hypothetical temperature observations given In table 1. 

The computed values of T^l are summarized In the lower half of 
table 1. 



given 



thus 

and 



To = 10 



= 2 

n = 4 



9 ; = i^)' = (- 2 ): > 

(-2); + 7i . 
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Table 1 



Some hypothetical temperature observations (Ti ) and 
the resulting Lamont-corrected temperatures (Tft ) 



Hour of observation (GOT) 00 06 12 18 24 

i 0 12 3 4 



Ti 

Si 

TiL 



2 3 6 9 10 



0 -2 -4 -6 -8 



2 12 3 2 



It must be noted that the temperatures as corrected by 
the Lament procedure are quite different than the original, 
uncorrected temperatures. However, since the present investi 
gation deals only with the meteorological temperature cycle 
(i.e. the T, -Lvalues), actual values of temperature are not im 
portant. 

If the actual temperature values were desired, it would 
be necessary to maintain a table of corrections, as well as 
the corrected values. 
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Theoretical discussion of the method used to 
separate the meteorological diurnal variation 
and the statistical variation 

As previously explained, the apparent diurnal temperature 
variation observed during the polar night for special types of 
days may be considered to be the sum of the meteorological var- 
iation and an effect resulting from statistical bias. More- 
over, the major purpose of this study is to determine an esti- 
mate of these two components. 

In order to illustrate how this dichotomy may be performed, 
it will be helpful to consider first how a particular observed 
dally temperature variation is "built up" by combining a 
meteorological influence with a statistical one. After the 
total diurnal variation has been constructed, it will once again 
be separated into its original components. This example will 
serve to Illustrate how the actual data for this study may be 
treated to derive the meteorological and statistical variations. 

First, consider that the set of T;l values listed in 
table 1 are representative of the meteorological variation. 

These temperatures will remain fixed for the particular hours 
for which they were reported (e.g., the meteorological influ- 
ence for 0600 GOT will always be IF, regardless of the starting 
time of the period) , 

Next, consider a purely hypothetical statistical vari- 
ation as listed in table 2. As previously pointed out by 
Hisdal, the statistical bias is not fixed at a particular 
value for a particular observation time, but depends upon the 
starting time of the period selected. To Illustrate this. 
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assume first that the starting time of the 24-hour period in 
table 2 Is 0000 GOT. Then the statistical bias for 0600 GOT 
would be IP, since it would be the second observation reported. 
On the other hand, if 0600 GOT were the starting hour of the 
24-hour period, the 0600 GOT bias would be 3P. Thus, the sta- 
tistical Influence may be thought of as a wave which is free 
to move through the diurnal or meteorological wave with con- 
tributions dependent upon the starting time of the 24-hour 
period. 

Table 2 

Hypothetical statistical bias 
for each 6-hourly observation 
time in a 24-hour period 

Observation time i T emn er atur e ( F ) 



1st 3 

2nd 1 

3rd 0 

4th 1 

5th 3 



Now, using the meteorological Influence from table 1 and 
the statistical influence from table 2, the "total" tempera- 
tures will be computed for various starting times, as shown 
in table 3. 

The values of total temperature for the various 24-hour 
periods are summarized in table 4. The smaller numbers adja- 
cent to the total temperature are the meteorological and 
statistical contributions to the total, as follows. 
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Table 3 



Total temperature variation resulting from the statistical 
and meteorological variations from tables 1 and 2 
respectively, for various starting times 



Starting 

Hour; 

Statistical: 

Meteorological: 

Total; 



Starting 

Hour: 

Statistical: 

Meteorological: 

Total: 



time: 0000 GOT 
00 06 12 18 

3 10 1 

2 12 3 

5 2 2 4 

time; 0600 GOT 
06 12 18 00 

3 10 1 

12 3 2 

4 3 3 3 



00 

3 

2 

5 

06 

3 
1 

4 



Starting 

Hour: 

Statistical: 

Meteorological: 

Total: 



time: 1200 GOT 
12 18 00 06 

3 10 1 
2 3 2 1 

5 4 2 2 



12 

3 

2 

5 



Starting time: 1800 GOT 
Hour: 18 00 06 

Statistical: 3 1 0 

Meteorological: 3 2 1 

Total; 631 



12 

1 

2 

3 



18 

3 

3 

6 
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total 



g 

tm; where t = 

s = statistical influence 
m = meteorological influence. 

The total temperature variations in table 4 are typical 
of those obtained from the GDC 1604 computer program. They 
are plotted in figure 3 in order to show the effect the stat- 
istical Influence has in shifting the time of the daily min- 
imum for various starting times of the hypothetical 24-hour 
period. 

Now, the next step in this discussion is to perform a 
dichotomy upon the total temperature values found in table 4 
in order to determine the original values of meteorological 
and statistical variation contributions. A more general proof 
of the separation method is given in section 5» Specifically, 
the method to be used is the same as that utilized by van der 
BI 3 I in his investigation of diurnal pressure variations. It 
consists of averaging the values of table 4 in either of two 
ways: by columns, to determine the meteorological influence, 
or diagonally, to determine the statistical Influence. Once 
either of these averages has been computed, the other may be 
determined by subtraction from the total values. These methods 
are shown below. The notation used is as follows: 

t = temperature (P) 

0 = m; refers to meteorological influence 

0 = s; refers to statistical Influence 

hh = hour of observation (GOT)* 
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Table 4 



Hypothetical total temperatures 
for various 24 hour periods 
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Figure 3: Hypothetical total temperature 
variation curves for various starting 
times of the period, based upon the 
values found in table 4. 
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Averaging by columns, and subtracting the constant value 
of the statistical bias, results in the meteorological varia- 
tion. 



^ LLLLOjJ.) 

4 /V 4 / 

= z 



4-^ 



= 1 



f L±J-±1±Q')- A - i. 

V 4 /■ 4 4 






= 13.^ 



= ^ 






IS 






= 3 



f / 3t I ->-0+ l \- tS _ S 






Averaging diagonally (assuming the 24-hour period begins 
at 0000 GOT), and subtracting the constant value of the mete- 
orological influence, results in the statistical influence. 

j. / £±JL1_£±-^ ) / 2.-«- i f 3. -*• 3 ) _ 20 g 

■^ 5.00 = ( 4 / \ ^ / 44 

= 3 



f / Z-^3 -h ^ t_2. ) .( I 3 *Jz) ^ IZ^ S. 

= 1 
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w =c- 
= o 



z+ 3 2. + I 







3 Z * I + Z 
4 




= 1 




2. •+ / + 2. •*- 3 

4 




= 3 



In this example we have recovered the values of the mete- 
orological and statistical Influences, assumed originally in 
constructing the values of total temperature. 

One difficulty in applying these methods to actual data 
lies in the fact that the constant terms which are subtracted 
from the averages of total temperature are not immediately 
apparent. This difficulty, however, is of no great concern, 
when the temperature trend and not the actual temperature per 
se is under Investigation. 

For the sake of completeness, the next section will 
illustrate how the use of these constants may be eliminated. 
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Elimination of the constant term, when computing 
the statistical bias 



The purpose of this section Is to Illustrate how the con 
stant term may be eliminated, when computing the statistical 
bias. The only known temperature values are those of total 
temperature, found In table 4, The hourly meteorological con 
trlbutlons are unknown. 

The notation used Is as specified In section 4, except 
for the following: 

Ojjj = a constant which Is dependent upon the 
meteorological Influence^ 

Cg = a constant which Is dependent upon the 
statistical bias. 

The column averages, computed In section 4, become: 
t = 12 - r 



Next, If the averaging Is done diagonally, using 0000 
GOT as the starting time, the following temperatures result: 




t - 12 - r 





tiYi,i-f r: Cs • 



5,06 ~ -f- , 



= ^ - Cm , 



0 . 



m 
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t - - r 

*'S,)8 ~ 4- '-''n J 



^ Cly, . 

Making the appropriate substitution for Cm in the last 
set of equations^ 

+■ _ ^ _ ^/Hjoo "tfn lA + 

*'5.00 — xf. 4 

£0 I / !1 * 3.4- 4- Ll\ ^ + 't.j o6 't 5yll tj 

= — 4 

f- _ £T — + ~t ijoo -t "tiyot + ts.iz + "65.18 . 

^.°o - /6 4 

- • ^^oo~^fy08~ ^S.u“'^iy/8 — 7 • (^) 

Slmilarlly, for the starting times 06, 12 and 18 h, we 
obtain, respectively: 

”^too“ ^iyi8 = 1 ^ (5) 



tsyOO -r tj-ofe ~36^/Z+ = 5 J 



(6) 



t$^oo+ isyO* 1 , (Y) 

Equations (4) through (7) comprise a set of four linear, 
homogeneous equations in four unknowns, which can be readily 
solved using the usual mathematical techniques. Such a solu- 
tion produces the following: 

tSyOO = 3 

"^S.Ofe = 1 

5.18 = ^ 

65^2.4 = tj^oo — 3 . 



19 



Thus, it is seen that if specific values of statistical 
bias or meteorological effect are desired, a set of homogeneous 
equations must be solved* However, as was stated in section 
4, this present study does not require specific temperature 
values, but rather their relative values* Therefore, the 
constants Cg and will not be removed from the computed 
diurnal averages* 
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General discussion of the actual results 



In general, the results of this study can be separated 
into four parts. The first is the set of 24-hour temperature 
variation curves, figures 4 through 25. These curves were 
drawn in order to determine to what degree statistical bias 
is apparent, when dealing with actual data. Also, in these 
curves the criteria used in selection of the type of day to 
be investigated were varied. The extreme cases considered 
were (1) wind speed greater than 15 kt, cloudiness greater 
than five-tenths; and (2) wind speed calm, cloudiness equal to 
zero. 

The second set of curves, figures 26 through 30 are the 
48-hour variations for various pairs of days selected. The 
selection criteria are comparable to those specified for the 
24-hour variation curves. The major purpose of the 48-hour 
curves is to show that the statistical bias effect is not res- 
tricted to a 24-hour period. 

The third set, figures 31 through 34_, comprise the 18-hour 
and the 72-hour variations. These curves are also included in 
the study to show how periods of various length effect the 
statistical bias. 

The fourth set of curves, figures 35 through 42, show the 
relationship between temperature variation and variations of 
atmospheric pressure, cloudiness, and wind direction. The pur- 
pose of this set is to aid in giving some explanation for the 
peculiarities found in the daily temperature trend. 

In all of these curves the abscissa is the hour of obser- 
vation in Greenwich Civil Time, which is the same time scale 
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that Hisdal used In his Investigation, It is convenient to 
use Greenwich Civil Time rather than local time, since the 
reporting times of surface observations are in GOT. 

If it is desired to obtain the diurnal variation, based 
on local time, add 11 hours to the observation times in Green- 
wich Civil Time, However, it must be remembered that local 
time loses its significance during the polar winter season, 
since the noon hour has no more sunshine than the "witching 
hour" , 
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The diurnal temperature variation for special 
selections of days 



Table 5 summarizes the comparison between the expected 
curvatures of the diurnal variations for selected types of 
day and various starting times of the period considered. The 
expected curvature is based on the discussion in section 1, 
There it was pointed out that days with high wind speed 
and/or cloudy skies will have above^normal average values of 
temperature, resulting in a convex diurnal temperature curve. 
Days with low wind speed and/or clear skies have below=normal 
mean temperatures, resulting in a concave diurnal temperature 
curve. In table 5 the shape of the daily meteorological 
course of temperature agrees favorably with the expected shape. 
However, as the criteria are made more restrictive, the agree- 
ment becomes worse, due to the increasing contributions 
arising from random effects as the sample size is reduced. 

(The effects of randomness are considered in section 11.) 

Table 6 contains a summary of the expected shape and 
actual shape of the statistical bias curves for the various 
types of day selected. The table also contains Information 
concerning the time of maxima or minima in the mean meteoro- 
logical curves and the range of temperature fluctuations 
involved. It is seen that the actual shape of the statistical 
bias curves agrees with the expected shape quite well. How- 
ever, as in table 5» making the selection criteria more 
restrictive produces more pronounced random errors, due to 
the reduction of the sample size. For instance, figure 21 
is based upon only seven days. The result is that no 
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Table 5 



Summarized results of the total temperature 
variation curves for various starting times 




3s Starting time = 1200 GOT 
4 s Starring time ss 1800 GOT 
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organized statistical influence is apparent in this case. In 
contrast, figure 7 is based on 174 days and the statistical 
bias has a definitely concave-shaped curve. 

One peculiarity that occurs in the set of temperature 
variation curves is the 2100 GOT maximum, appearing in the 
statistical bias curves. This maximum, which would occur at 
0800 local time, gives an unnatural appearance to the bias 
curves. Therefore, considerable time was spent in determining 
the cause of the peculiarity. One resulting conclusion is 
that the 2100 GOT maximum is associated with the criterion of 
calm winds, since it appears in all curves in which calm winds 
are included, except those of figure 21 , Here the maximum 
occurs at 1800 GOT, However, this latter figure is based on 
only seven days and one should be careful not to attribute 
undue significance to such a small sample. 

An attempt was made to exclude the calm days in computing 
the mean diurnal temperature variation in order to determine 
whether the 2100 GOT maximum would disappear from the statisti- 
cal bias. The results are contained in figures 5 and 23 » in 
which the maximum does not appear. However, the criterion for 
figure 23 (wind speed between two and 15 ht) was too restrictive 
and only four days could be used. 

Further investigation revealed how the calm-day observa- 
tions affected the statistical bias curves. This matter is 
explained in more detail in section 11 . It will suffice the 
present discussion to say that the maximum in the statistical 
bias curve resulted from calculating the bias through the use 
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of a total temperature curve which had an above-normal random 
contribution at 2100 GOT, 

It is Interesting to note that the 1 8-hour statistical 
bias curve, figure 32, does not contain a maximum comparable 
the the 2100 GOT maximum. Therefore, it is doubtful that such 
a maximum would occur with an independent selection of data. 
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Figure 4; Total temperature variation curves 
with variation in sWrting time of the period, 
(wind speed -^15 ht; cloudiness -is* 5/l0) 
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Figure 5: Comparison of the total temperature 

variation with the meteorological (Met.) and 
the statistical bias variations. 

(wind speed ^15 ht; cloudiness 5/lO) 
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Figure 6: Total temperature variation curves 

with variation in starting time of the period, 
(wind speed kt; cloudiness'^ 5/l 0) 
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Figure 1 : Comparison of the total temperature 

variation with the meteorological (Met.) and 
the statistical bias variations. 

(wind speed 1^15 ht; cloudiness 5 /l 0 ) 
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Figure 8: Total temperature variation curves 

with variation in starting time of the period, 
(wind speed 8 let; cloudiness 5/l0) 
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Figure 9- ComparisorL of the total temperature 
variation with the meteorological and the 
statistical bias variations, 

(wind speed ^8 kt; cloudiness^5/lO) 




-33 



T EMPERATURE (F; 



Figure 10: Total temperature variation curves 

with variation in starting time of the period, 
(wind speed ^5 kt; cloudiness^5/lO) 
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Figure 11: Comparison of the total temperature 

variation with the meteorological and the 
statistical bias variations. 

(wind speed ^5 ht; cloudiness^^S/^ 0) 
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Figure 12: Total temperature variation curves 

with variation in starting time of the period, 
(wind speed ^4 kt; cloudiness 5/l0) 
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Pigure 13: Comparison of the total temperature 
variation with the meteorological and statis- 
tical variations. 



(wind speed -4 kt; cloudiness;^5/lO) 
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Figure 14: Total temperature variation curves 

with variation in starting time of the period, 
(wind speed --3 kt; cloudinessi5/l0) 
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Figure 15: Coinparison of the total temperature 
variation with the meteorological and 
statistical bias variations. 



(wind speed 3 ht; cloudiness 5/1 0) 
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Figure l6: Total temperature variation curves 

with variation in starting time of the period, 
(wind speed 2 kt; cloudiness^5/l 0 ) 
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Figure 17: Comparison of the total temperature 

variation with the meteorological and the 
statistical bias variations, 

(wind speed 2 kt; cloudiness^5/lO) 
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Figure 18; Total 
with variation in 
(wind speed = 0 



temperature variation curves 
starting time of the period, 
kt; cloudiness 5/l 0) 
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Figure 19: Comparison of the total temperature 

variation vdth the meteorological and the 
statistical bias variations. 

(wind speed = 0 kt; cloudiness=5/l0) 
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Figure 20: Total temperature variation curves 

v.'ith. variation in starting time of the period, 
(wind speed = 0 ht; cloudiness = 0} 
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Figure 21: Comparison of the total temperature 

variation with the meteorological and the 
statistical bias variations, 

(wind speed = 0 kt; cloudiness =0) 
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Figure 22: Total temperature variation curves 

with variation in starting time of the period. 
(2 kt-i. wind speed ^15 kt; cloudiness^ 5/l0) 
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Pigure 25' Comparison of the total temperature 
variation with the meteorological and the 
statistical bias variations. 

(2 kt — wind speedih.15 kt; cloudiness rS;. 5/lO) 
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Figure 24: Total temperature variation curves 

vrith variation in starting time of the period, 
(wind speed U-.15 ht; cloudiness r^ilO/'iO) 
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Pigure 25: Comparison of the total temperature 
variation with the meteorological and statis- 
tical variations, 

(wind speedh-15 ht; cloudiness A- 1 0/l 0) 
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8. 48-hour temperature variation 

As previously mentioned » the purpose of the 48-hour tem- 
perature variation analysis is to examine whether the statis- 
tical bias occurs for periods other than 24 hours. The results 
of this analysis are presented graphically in figures 26 
through 30. 

Figure 27 clearly shows an effect of statistical bias. For 
this figure the selection criteria were windy and/or cloudy 
periods, which, according to Hisdal, should have their highest 
temperatures near the middle of the period. This is what ac- 
tually occurred, but only five 48-hour periods met the speci- 
fied selection criteria. 

In figure 29, a total of 39 periods met the criteria, but 
the actual results are inconsistent with those expected. Accor- 
ding to the discussion of section 1 , for days with wind speed 
less than or equal to 15 knots and cloudiness less than or 
equal to five-tenths the temperature trend should be concave 
in shape. In figure 29? the statistical bias is convex. This 
departure from a priori expectation may be due to the random 
variations included in the total temperature curve, from which 
the bias was deduced. This idea is treated more thoroughly in 
section 11, 

The last figure in the 48-hour set, figure 30, Illustrates 
the mean 48-hour temperature variation which results when no 
particular selection criteria are specified. In theory these 
curves should have negligible statistical bias, considering the 
discussion in section 1 . This theory is essentially borne out 
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in figure 30, where the mean meteorological variation, obtained 
by the column averaging method. Is very nearly the same as the 
total variation. Thus, upon subtraction of the meteorological 
from the total variation, the resulting statistical bias does 
not appear significantly different from zero. 

It should be noted that the conclusion of a zero bias 
through such a procedure is complicated by end effects due to 
different starting times. For example, the first two three- 
hourly temperatures used In computing the 0000 GOT to 0000 GOT 
mean curve are not included In the 0600 GOT to 0600 GOT curve. 
Likewise, the last two temperatures In the 0600 GOT to 0600 GOT 
curve are not Included In the 0000 GOT to 0000 GOT curve. 
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Figure 26: 48-hour temperature variation 
curves with variation in starting time 
of the period. 

(wind speed ^15 ht; cloudiness ^5/lO) 
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Figure 27 J Comparison of the total 48-hour 
temperature variation with the meteorological 
and statistical bias variations. 

(wind speed =^15 ht; cloudiness^ 5/10) 
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Figure 28; 48 -tLOur temperature variation curves 
with variation in starting time of the period, 
(wind speeder 15 kt; cloudiness:^; 5/l0) 




wQur.vei 



'-itJhfttfl 



-0u^v^2^ 



iTTffica se 



Legend: 

Basic curve: — 
reoeated curve 









■0urve4' 



00 06 12 18 00 06 12 I8 00 06 

HOUR OF OBSERVATION 
(GOT) 



nir. 

ni r 


pd-t 

4|tl 


ijil: 


iii 



- 54 



temperature: (f; 



Figure 29 : Comparison of the total 48-hour 
temperature variation with the meteorological 
and statistical bias variation. 

(wind speed^15 ht; cloudiness <=5/10) 
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Figure 30i 48-Hour temperature variation 
curves with variation in starting time 
of period(for all observations, regard- 
less of wlndspeed and cloud amount). 
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9. 1 8-hour and 72-hour temperature variations 

As with the 48-hour temperature variation, the purpose 
for studying the 18-hour and 72-hour variations is to demon- 
strate that a statistical bias effect is not limited to a 
24-hour period. 

The 1 8-hour temperature variation curves, figures 31 and 
32, clearly depict the existence of a statistical bias. In 
principle, these curves should be concave, since the criteria 
specified are wind speed less than or equal to 15 knots and 
cloudiness less than or equal to five-tenths. The actual 
curves are, in fact, concave. 

In figure 32 the meteorological temperature variation and 
the statistical bias have been computed by means of the two 
schemes proposed in section 4; that is, by both the column- 
averaging and diagonal-averaging processes. The resulting 
curves are quite similar. 

Figures 33 and 34 depict the 72-hour temperature varia- 
tion curves. These curves are designed to test the hypothesis 
that the statistical bias effect exists for considerably 
longer periods than 24 hours. However,, they do not support 
this hypothesis. In fact, there is evidence of an upper limit 
for the length of the period yielding a bias effect. A possi- 
ble explanation for this limitation is that the physical 
mechanism causing the type of day to persist for a long period 
remains long enough so that the observations on either side of 
the selected period also meet the criteria. If this were the 
case, no statistical bias would be evidenced. 



57 







I 

I 



i 

I 

I 




A further investigation is necessary in order to deter- 
mine which synoptic situations at Antarctica are associated 
with long sequences of similar days. 

One difficulty which arises in studying a period as long 
as 72 hours is that the number of periods meeting the criteria 
are relatively few. Thus, the statistical significance is 
greatly reduced compared to that of shorter periods. 
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Figure 31 1 1 8-hour temperature variation 
curves with variation in starting time 
of the period. 



(wind speed:Sl5 kt; cloudinesss55/lO) 
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Figure 32: Comparison of c,he total 1 8-hour 
temperature variation with the meteorological 
and statistical bias variations. 

(wind speed -615 kt; cloudiness 2 ^ 5 / 10 ) 
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gure53: The 72-hour temperature variation with variation in starting 
me of the period, (wind speed^15 kt; cloudlness:^/lO) 
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Figure 34; Comparison of the total 72-hour 
temperature variation with the meteorological 
and statistical bias variations, 

(wind speedi^ 15 ht; cloudiness^ 5/10) 
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10. The diurnal temperature variation, based upon 

all observations 1956 - 1961 

Figure 35 Is a chart of mean diurnal or meteorological 
temperature variation curves, based on all three-hourly sur- 
face observations at McMurdo Sound, during the dark season, 
for the years 1956 through 1961. The first six of these curves 
are annual means, based on observations of 93 days. The last 
curve is the six-year mean diurnal variation, based upon 558 
days. 

Since no selection criteria were specified for these mean 
curves, each represents a diurnal temperature variation, that 
is, one free from statistical bias. 

There are two pronounced features of the six-year mean 
curve. The first is the temperature maximum occurring at 1500 
GOT. This maximum, however, may have no particular significance, 
since it is caused, in part, by the abnormally high 1500 GOT 
mean temperature which occurred in 1956. Moreover, during the 
years 1958 through 1961, there is a minimum value at 1500 GOT. 

The second feature of the six-year mean diurnal temperature 
variation curve is the maximum occurring at 0000 GOT. This 
appears to be a regularly occurring phenomenon, at least for 
the years included in this study. 

In order to discover a possible explanation for the 0000 
GOT maximum, the diurnal variations of some other meteorolo- 
gical parameters were plotted in figures 36 through 42. The 
first of these, figure 36, gives an estimate of the "mean diur- 
nal variation of wind direction," In order to simplify the 
computation, a scalar-average wind direction was computed 
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Instead of a vector average* However, the diurnal course of 
wind direction did not reveal a modal value at 0000 GOT* On 
the other hand, Hisdal [1, p,172] in his Investigation did, 
in fact, obtain a relationship between wind direction and sur- 
face temperature, using a histogram of wind direction in con- 
junction with the diurnal temperature trend. It would be 
interesting to consider histograms of the wind-direction fluc- 
tuation associated with the temperature curves of this study. 

The next parameter examined was the diurnal variation of 
cloudiness. The annual means and six-yearly mean are found 
in figure 37. This set of curves bears a very close resem- 
blance to the mean curves of diurnal temperature variation, 
even to the point of having a 0000 GOT maximim. Once again, 
the conclusion may be drawn that cloudiness plays a very im- 
portant role in the mean diurnal temperature fluctuation during 
the polar night. However, it cannot be concluded that the 
diurnal fluctuation of cloudiness is the cause of the diurnal 
temperature variation, for then it becomes quite difficult to 
explain why there should be a systematic periodicity in cloud- 
iness. It seems more reasonable to belie\’’e that the two 
phenomena are due to some common cause. 

A relationship between cloudiness and temperature would 
also be indicated, if a comparison were made between figures 39 
and 19. In these curves, the selection criteria were calm wind 
and cloudiness less than or equal to five-tenths. These cri- 
teria were chosen in order to minimize temperature changes 
due to turbulent mixing associated with higher wind speeds. 
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Comparing the mean variations in figures 19 and 39 » it appears 
that for periods with calm winds a maximum or minimum of tem- 
perature tends to lag the corresponding cloudiness by three 
hours. Table 7 summarizes some observations of this lag re- 
lationship. 

Table 7 

Time (GOT) of occurrence of maxima/minima 
for days with wind speed = 0 and cloudiness 
less than or equal to five-tenths 

Cloudiness Temperature 



Max, 


Min. 


Max, 


Min, 


2400 




2100 






2100 




1800 


1800 




1500 






1500 




1200 


1200 




0900 






0900 




0600 


0300 




0300 




neither 


mean wind 


direction nor 


mean cloudiness 



yielded a causative explanation of the diurnal variation of 
temperature, another parameter was investigated. The parameter 
selected next was the sea-level pressure. Its annual mean 
diurnal variation and six-year mean are found in figure 40, 
Also, figure 42 shows a comparison of the six-year mean diur- 
nal pressure variation with that of the corresponding temper- 
ature variation, as well as the variation of other parameters. 
The Important features of the slX”=year mean diurnal 
pressure variation is that it has a maximum value at 1200 GOT 
and a secondary maximum at 0000 GCT, This result agrees with 
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the results of a recent Investigation of Antarctic data by 
Carpenter [10] , He used two years of three-hourly surface 
pressure data In making a harmonic analysis to determine the 
phase and amplitude of the diurnal and semidiurnal pressure 
waves. Carpenter’s results for McMurdo Sound are listed in 
table 8 which compares favorably with the pressure variation 
In figure 42. 

Table 8 

Mean harmonic components of the diurnal 
pressure wave at McMurdo Sound, after 
Carpenter [10] 

Ampl.(mb) Mean time of max. 
Diurnal Oscillation .001 2118 GOT 

Semidiurnal Oscillation .076 1040 GOT 

In view of the fact that there is a definite diurnal and 
semidiurnal pressure fluctuation at McMurdo Sound, the next 
problem is to determine what Influence they have, if any, on 
the diurnal temperature variation. 

Since the semidiurnal variation has its maximum value 
near 1200 GOT, it is readily apparent that this compressional 
pressure wave can not account for a temperature minimum which 
occurs at the same time. The Poisson equation shows that for 
an adiabatic pressure increase, a temperature increase will 
result. Moreover, the diurnal pressure wave (maximum at 2118 
GOT) has an amplitude of only .001 mb. This pressure increase 
would result in a temperature increase by adiabatic compression 
of the order ,0002P, 

Thus, it is seen that neither the diurnal nor semidiurnal 
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?igure35: Mean diurnal temperature 
variation during the polar night 
for the years 1956 through 1961. 
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V.’IMD DIRECTION 



Pigure3o: Kean diurnal variation of vrind 
direction during the polar night for the 
years 1956 through 1961, based upon all 
observations, regardless of wind speed 
and cloud amount. 
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Pigure 37 : Mean diurnal variation of cloudiness 
during the polar night for the years 1956 
through 1961 , based upon all observations, 
regardless of wind speed and cloud amount. 
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CLOUDINESS (tenths) 



Pigure 38: Total cloudiness variation 
curves with variation in starting time 
of the period. 

('vn.nd speed = 0 kt; cloudiness:^ 5/l0) 
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Ptgure 39 ‘ Comparison of the total 
cloudiness variation with the meteor- 
ological and statistical "bias variations, 

(wind speed = 0 let; cloudiness:^: 5 / 1 0) 
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?igure40: Mean diurnal variation of atmospheric 
pressure during the polar night for the years 
1956 through 1961 5 based upon all observations, 
regardless of wind speed and cloud amount. 
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Figure 41: Comparison of the total 
atmospheric pressure variation with 
the meteorological and statistical 
bias variations. 



(wind speed = 0 kt; cloudiness = 0) 
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Figure 42: Diurnal variation of temperature, 
atraospheric pressure, wind direction, and 
cloudiness during the polar night. These 
curves are "based upon six winters of 
o"bservations, 1956 through 196 I. 
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pressure variations can account for the 0000 GOT maximum in 
the diurnal temperature variation. The diurnal pressure wave 
is not of sufficient magnitude to affect the diurnal tempera- 
ture variation, and the semidiurnal pressure wave is 12 hours 
out of phase with the temperature. 

With the foregoing physical factors eliminated as rea- 
sonable explanations for the 0000 GOT maximum in the daily 
temperature curves considered previously, it is necessary to 
again consider a statistical explanation, 

A climatological statistical effect has been proposed as 
an explanation by van der Bljl^who suggested that the seasonal 
temperature trend could produce such a maximum. A schematic 
version of the seasonal effect is shown in figure 43. 



Figure 43 

Schematic Illustration of 
the seasonal trend effect 
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A 24-hour segment of the seasonal temperature trend, 
denoted AB, is enlarged and shown as a dashed line in the 
upper right-hand corner of figure 43» The solid lines in the 
enlarged portion show the segment as corrected by the Lament 
procedure. 

If the entire temperature trend for the cooling season 
were divided into 24-hour segments, each would have a concave 
curvature, because the over-all trend is concave. Now, if 
these 24-hour segments were averaged for the entire cooling 
season, the mean would also be concave. This concavity would 
give the appearance of a diurnal temperature fluctuation, such 
as that observed in the present investigation. • 

Figure 35 shows that the concavity of the mean diurnal or 
meteorological temperature variation produces a temperature 
range of O.SF. It is of interest then to determine whether 
the seasonal temperature trend actually has enough curvature 
to produce 0.5F in a 24-hour period. In pursuing this inves- 
tigation, the seasonal temperature trends were plotted for 
each of the years 1956 through 1961. Also, the six-year mean 
trend was computed. These curves are shown as figures 44 
through 50. 

Examination of the seasonal temperature trends reveals 
that they are not as smooth as was expected. Also, the cur- 
vature from the beginning to the end of the dark season is 
not of sufficient magnitude to cause a five-tenths degree 
temperature range in a 24-hour segment. Therefore, a temper- 
ature-variation model which would increase the curvature of 
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the average 24-hour segments would substantiate van der Bljl's 
proposal and yet produce, upon averaging, the long-term mean 
seasonal temperature curve. 

This curvature Increase may be accomplished by considering 
that actual observations of the seasonal temperature trend In a 
given year are as shown In figure 51. In this figure the 
dashed line represents the long-term seasonal trend, which Is 
the same as In figure 43. The solid line represents the 
unsmoothed temperature trend, which shows that relatively long 
cooling periods during the dark season are Interrupted by short 
periods of above seasonal temperature. 

The effect Is to produce a series of “spiked" temperature 
maxima which greatly Increase the mean curvature of the 
average 24-hour segments during the prolonged cooling segments. 
With this refinement the mean diurnal temperature variation 
may have enough curvature to produce a range of 0.5F. 

A more complete study of the seasonal temperature trend 
Is necessary, however, before Its effect on the curvature can 
be validated. In particular, If the shape of the summer tem- 
perature trend Is as shown In figure 51 , and the total mean 
diurnal temperature variation for the mid-summer months has a 
minimum at 0000 GOT, then the seasonal trend effect would have 
some merit for the summer months also. 
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Figure'^H: Seasonal trend of the temperature during the polar 
night at McMurdo Sound, Antarctica. (15 May - 15 August, 1956) 
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Figure "^-5; 'Seasonal trend of the temperature during the polar 
night at McMurdo Sound, Antarctica. (15 May - 15 August, 1957) 
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Figure Seasonal trend of the temperature during the polar 
night at McMurdo Sound, Antarctica, (15 May ~ 15 August, 1958) 
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Pigure^f/': Seasonal trend of the temperature during the polar 
night at McMurdo Sound, Antarctica. (15 May - 15 August, 1959) 
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Figure 50; Mean seasonal trend of the temperature during the 
polar night at McMurdo Sound, Antarctica. (1956 - I961) 




o 

4- 



IC\ 

+ 



LA 

1 



O 

8 

id) 



a 

C\J 

s 



Lf\ 

I 

STdfllVHEdKai 



A 

OJ 

I 



o 

A 

8 



A 

A 

I 



O 

I 



A 

I 



84 



Figure 51 



Schematic model of a proposed seasonal 
temperature trend for the polar regions 
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1 1 , Discussion of random errors 

Throughout this study, random errors have been mentioned 
as an Important factor In the erratic behavior of various 
curves. For Instance, the 2100 GOT maximum In the statistical 
bias curve of figure 7 was considered to be caused by the 
randomness of some of the data used In computing the mean 
values. The following discussion will Illustrate that the 
2100 GOT was. In fact, caused by random errors. Also, a com- 
putational scheme will be suggested, which reduces the effect 
of these errors. 

The notation used Is as follows: 

M] = temperature due to the meteorological 

Influence at observation time 1 = 0,,,, n,,,,4 
Sj = temperature due to the statistical bias 
t„^ = total temperature = + Sj + 

m = 0, 4 

n — 0,,,,,3 

e^,,= temperature due to the random error 

cr^ = variance of the single- observation random 

2 - 

error = . 

Subscripts: 

m = row In matrix 
n = column In matrix , 

1,3= chronological order of the temperature 
observation . 

Consider table 9 to be a matrix of total temperature 
similar to table 3. The meteorological effect Is constant 
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for any particular column, and the statistical bias is con- 
stant for any particular diagonal row. The errors are assumed 

/ 

to be random and uncorrelated in both the vertical and the 
diagonal directions. 



Table 9 

Generalized Theoretical Table 
of Mean Total Temperatures 




Using the column-averaging method, as explained in sec- 
tion 4, the following totals result from using the temperatures 
in table 9- ^ 

Column 1 : 4Mn + 

rn-o 



Column 2: 4 ^ S; + ^ 

J'O 

Column 3: 

Column 4: 4 ^^3 * 

J'O 

Next, divide by four to obtain the column mean temper- 
atures: 

,3 3 

Column 1; No + 4 X 4 X 

j=^o 

Column 2: N, + { ^ 5j + ^ X^^i 

j-O 
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Column 3 5 
Column 4: 






' ^ro 



•>M = 0 



Subtracting the column averages from their corresponding 
first-row entry yields an estimated value of the statistical 
bias for the various columns; 



Column 



Column 



Column 



Column 



1 : 


(M. 
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s. + 


e.„) 

hi:0 
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- i (S, + S 


S 3) + ^ 600 




^zo ^30) 


2: 


(M, 


+ S, -b ( 
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-(M, 




' m-o / 








X Q 

1 








te^, -nej) 


3 : 


(M. 


+ Sj. + e 


ox ) ^ 








~(Mz 


+ i 


4 4 ^ ^‘^'0 








^ s 

■t 


-i(s,^S 


. + S3) t ^ 


0 - 


12 


4 ; 


(Mj 


+ S3 + e 


0?) 








• (m. 




r ^ ^ nmj 

^ m*o / 







I + S, ^Si)+ I e ,3 + 



The next step will be to determine the variance of the 
statistical bias, in order to be able to compare it with the 
variance of the statistical bias computed by diagonal-aver- 
aging. 

If o’ is the variance of a single-observation random error, 
then the variance of the first column- average equals 
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The variance of the statistical bias is then computed as 



are considered to be fixed numbers and therefore have zero 
variance. 



The other columns treated similarly will yield the same re- 
sult, which is the standard deviation of the 

The second method considered in section 4 is that of 
averaging the total temperatures diagonally. As explained 
previously, a diagonal average will produce the statistical 
bias, without having to subtract the average from the total 
temperature. As will be proven below, this method results 
in a more accurate estimate of the statistical bias than that 
arrived at through a column average. Using general notation 
we have 

Diagonal sum; 4 Sj X ^ ^ ; j 



3 



shown below. Note in this connection that Mq and 




For column 1 : Var [ ( ^>°)j 

~ ^ V(^r too ^ 7t (’Vftirejo+l/c/rejo + V^r e,p) 




Average; 
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Variance in average; 



1/ Var ^ 

i-i 

- ]6 X til 

- h 

= ir". 

Note that the diagonal method reduces the variance in the 
statistical bias by a factor of three » compared with the var- 
iance resulting from the column method. . 

An important point that must be emphasized is that the 
reduction factor of three is valid only when the temperature 
observations are six hours apart. In other cases, where the 
observations are other than six hours apart, different re- 
duction factors are v'alid, 

A generalized derivation, similar to the one for the 
six-hour spacing, can be used to arrive at the values in 
table 10, which are the factors by which the variance in the 
statistical bias will be reduced if a diagonal average is used 
in computing the bias, rather than a column average. 

Table 10 

Reduction factor, R, for the error 
variance in the statistical bias 

At 1 2 5 4 6 12 24/k 

R 23 11 7 5 3 1 k-1 

At = time between observations 

- (24/k) hours 

R = Reduction factor 

k — 1s2,3,oo.o9n 

This discussion of random errors has shown that the best 
method for obtaining the statistical bias is by averaging the 
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temperatures In table 9 diagonalljo A similar discussion 
would also prove that the best method for obtaining the meteor- 
ological influence is through a column average. 

Figure 52 was plotted in order to provide a graphical 

/ 

comparison of the results of these two methods. The most 
striking feature of the figure is that using the diagonal 
average in computing the statistical bias has eliminated the 
2100 GOT maximum. 
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Figure 52: Comparison of the total temperature 
variation v:ith the meteorological and statis- 
tical variations. 

(vrind speed lA 15 ht; cloudiness 5/lO) 




09 12 15 18 

HOUR OF OBSERVATION 
(GOT) 



92 



I 

■1 



12. Conclusions and recommendations 



The purpose of this Investigation was to determine the 
nature of the diurnal temperature variation during the polar 
night at McMurdo Sound, Antarctica. The conclusion based upon 
data for the years 1956 through 1961 Is that no diurnal tem- 
perature variation, In the normal sense, exists during the 
polar night. The maxima that appear In the dally temperature 
trend are due to synoptic fluctuations. The concave curvature 
between temperature peaks appears to reflect that of the sea- 
sonal cooling curve. 

It Is recommended that the column-averaging procedure 
should be used only for determining the meteorological compon- 
ent of the diurnal temperature variation. For computing the 
statistical bias, the diagonal-averaging method should be used. 
With this scheme, the variance In the statistical bias can be 
minimized. 

For further study It Is suggested that Investigations, 
similar to the present one, be made for the following situations 

(1) summer season, when the sun never goes below 
the horizon^ 

(2) other polar stations, both In the winter and 
summer seasons. 
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